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' We present a detailed investigation of the particle pair separation process in homogeneous isotropic 

, turbulence. We use data from direct numerical simulations up to Rx ~ 280 following the evolution 

04 ■ of about two million passive tracers advected by the flow over a time span of about three decades. 

We present data for both the separation distance and the relative velocity statistics. Statistics 
I are measured along the particle pair trajectories both as a function of time and as a function of 

their separation, i.e. at fixed scales. We compare and contrast both sets of statistics in order to 
gain an insight into the mechanisms governing the separation process. We find very high levels of 
intermittency in the early stages, that is, for travel times up to order ten Kolmogorov time scales. 
The fixed scale statistics allow us to quantify anomalous corrections to Richardson diffusion in the 
inertial range of scales for those pairs that separate rapidly. It also allows a quantitative analysis of 
' intermittency corrections for the relative velocity statistics. 



U 



o 



PACS numbers: 47.27.-i 47.10.+g 



INTRODUCTION 



> 

, The relative dispersion of pairs of particles is important because of its connection with the problem 
^) ' of concentration fluctuations 0,0,13 and because of the insight it provides into the spatial structure 
of turbulent flows. In contrast with single particle dispersion, which is mostly driven by the 
' large scale - energy containing - eddies, the dispersion of pairs of particles depends on velocity 
fluctuations of order the separation of the pairs. Thus, the early stages of relative dispersion, up 
to the integral time scale, are expected to reflect the universal nature of small scale turbulence 
^ (independent of the large scale flow) and the intermittent character of the energy cascade. The 
. latter appears to manifest itself in some particle pairs remaining close together for long periods of 
S ' time while others separate rapidly. 

> 
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Clearly, a good understanding of the mechanisms of relative dispersion will lead to better models. 
Among key features of relative dispersion are a separation dependent time scale and long-time 
correlations of quantities such as the relative velocity. Many different types of quantitative models 
of relative dispersion have been proposed including Lagrangian stochastic models e.g. |^^I3 and 
kinematic simulation . For a review of relative dispersion and Lagrangian stochastic models 
in particular, we refer the reader to Sawford In recent years it has become apparent that the 
success of these models at small scales will depend on their ability to capture the intermittency 
of the separation process |^. One purpose of this paper is to provide a detailed quantitative 
and qualitative analysis of the separation process which we hope will eventually lead to improved 
models. 



Results from observations of the spread of marked particles (pairs or clouds of tracers in the atmo- 
sphere and in the ocean), summarised in classical textbooks such as those by Monin and Yaglom 
and Pasquill and Smith p^ . is testimony to the difficulty in getting reliable experimental data in 
fully developed turbulence. Although much prog ress has been made in recent years in experimental 
measurements of single Lagrangian particles |llLll2l |. relatively few Lagrangian measurements have 
been obtained following pairs of particles. A notable exception is Ott and Mann p^'s] who report 
Lagrangian inertial range scaling even at modest Reynolds numbers, of the order R\ ~ 100, where 
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i?A is the Taylor scale Reynolds number. As a result, direct numerical simulation (DNS) is still 
the most important source of detailed Lagrangian statistics (e.g. 0, 0, 0, 0, at Reynolds 
numbers up to order R\ ~ 280). In this paper we analyse the results of a recent DNS of 3 — D 
homogeneous isotropic turbulence seeded with Lagrangian particles '19] . Although homogeneous 
isotropic turbulence has limited application to real situations, it is the simplest configuration for 
studying the statistics of relative dispersion. 

The paper is organised as follows. In section^ we outline the numerical scheme used for calculating 
the data and discuss statistical uncertainty and variability. We present results on the statistical 
properties of both the particle pair separation and its relative velocity. These are considered in 
sections IlIII and HVI respectively. In both cases we compute the statistics as a function of time and 
as a function of separation, that is, at fixed scales. The latter allows for a more accurate separation 
of the dissipative, inertial and the integral scale regions. 



II. DNS METHODOLOGY 



The direct numerical simulation of homogeneous isotropic turbulence was performed on 512'^ and 
1024'^ cubic lattices with Reynolds numbers R\ ^ 180 and R\ ~ 280, respectively. The Navier- 
Stokes equations were integrated using fully de-aliased pseudo-spectral methods for a total time 
spanning nearly three decades (from the order of a tenth of the Kolmogorov time scale, r^, to 
approximately three times the integral time scale, Tt,). The flow was forced by keeping the total 
energy constant in the first two wavenumber shells [23|. The flow at Rx — 284 was seeded with 
approximately two million Lagrangian passive tracers once a statistically stationary velocity fleld 
had been obtained. The positions and velocities of the particles were stored at a sampling rate 
of 0.07r^. The numerical parameters are summarised in Tabled In this DNS, dissipative scales 
are well resolved, satisfying rj ~ dx, where da; is the grid spacing. The Lagrangian velocity was 
calculated using linear interpolation which was demonstrated to be adequate for obtaining well 
resolved particle accelerations pll |. 

The particles were initially arranged in tetrads which were uniformly distributed in the flow. A 
total of 960,000 pairs with initial separations tq = 1.2i] and ro = 2.577 were formed this way. 
Particle pairs with larger initial separations were formed by following pairs chosen from different 
tetrads. In this way we also follow pairs with initial separations tq = 9.8?7 and rp = 19.6?]. The 
number of pairs varied from 5.10^ to 1.10^ depending on the chosen initial separation. In both 
cases a particle may have been used more than once to form a pair. 

Previous studies have shown that Lagrangian statistics are affected by highly non-Gaussian fluc- 
tuations (see e.g. 0,^3)- Thus, it is important to quantify the statistical uncertainty of some 
typical variable in order to ensure the reliability of the results within our statistical sample. Sta- 
tistical errors were estimated by dividing the sample into five sub-ensembles and calculating the 
minimum and maximum values. We find that the error is at worst approximately 15% for the 
separation skewness and approximately 25% for the relative velocity skewness. The nature of the 
forcing scheme used in the present DNS meant that relatively little temporal variability of globally- 
averaged quantities was observed (see Overholt and Pope for discussion on forcing schemes 
and temporal variability). In particular, fluctuations about the mean of the energy dissipation, e, 
were at most 10% during the evolution of the DNS. Thus, in the following we may safely use £ 
(and other globally-averaged quantities) to scale the two-particle statistics. 
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III. SEPARATION STATISTICS 



A. Fixed time statistics 



We consider the motion of two marked fluid particles, labelled by the superscripts (1) and (2). In 
homogeneous turbulence, it is sufficient to consider the statistics of the instantaneous separation 
of the positions of the two particles, namely r{t) — r^^\t) — r^'^\t). Furthermore, in isotropic 
turbulence, the separation magnitude r = |r| plays a fundamental role in the problem of relative 
dispersion. 

Following the well known ideas of Richardson j23l | , relative dispersion in the inertial range of time 
scales, <^ t ^ Tl, can be described in terms of a diffusion equation for the probability density 
function (pdf) of the pair separation p(r, t). In spherical coordinates this is given by 

dt r2 dr V dr J ' ^ ' 

where K{r) is a scalar eddy diffusivity. On the basis of experimental measurements in the at- 
mosphere, Richardson proposed that K{r) = kos^^^r'^^^ where fco is a dimensionless constant. 
Assuming a small enough initial separation and a large enough travel time, it can be shown (see 
e.g. Monin and Yaglom |^, p. 574) that a spherically symmetric solution of is given by 

where A = (3/2)^/r(9/2) is a normalisation constant. This exhibits strong non-Gaussianity with 
a narrow peak at the origin and very large tails and gives rise to the celebrated scaling for the 
second order moment 

(r^) = get\ (3) 

Here g — 1144fco/81 is the Richardson constant which is supposed to be universal. This result was 
also derived by Obukhov [2^ using Kolmogorov's classical theory of turbulence (K41) 



The Richardson pdf is perfectly self-similar, all positive moments behave according to the di- 
mensional law: rP oc t^^/'^. The scaling Q is notoriously difficult to achieve both in laboratory 
experiments and in DNS on account of the large separation of scales that is required to observe it. 
As a result, estimates of g have varied widely, from 0.06 to 3.5 8]. The main practical difficulties 
in achieving a long inertial subrange are due to dissipative range effects at the ultraviolet end of 
the spectrum and integral scale effects at the infrared end of the spectrum. In the dissipation 
range, pairs separate exponentially and with widely varying growth rates - some pairs separate 
rapidly while others remain close together. This leads to the formation of a broad distribution 
of separations. As a result, slowly separating pairs (which remain in the dissipative range) and 
rapidly separating pairs (which approach the integral scales) 'contaminate' the statistics in the 
inertial range. A very large Reynolds number is therefore required to produce reliable Lagrangian 
statistics in the inertial range. 



In Figure n we plot the mean square separation (r^) versus normalised by the Kolmogorov 
microscales, r] and respectively. Although the curves begin to collapse at large t, they do not 
display a scaling and still show a dependence on the initial separation. Thus, any attempt to 
extract the value of the Richardson constant will be marred by the memory of the initial separation. 
The simplest way to measure g is to plot (r^) scaled by the asymptotic prediction, et^ , and look 
for a plateau. These curves are displayed in the inset of Figure ^ It is clear that none of them 
produces a good plateau and given the spread of curves with different initial separations, the value 
will be at best an order of magnitude estimate subject to considerable uncertainty. 

An alternative method, used in 0,0,0], consists of fitting a straight line to (r^)^/"^ in a suitable 
time interval. If equation Q holds, this straight line, when extrapolated back towards t = 0, 
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should pass through the origin and have a slope of (ge)^/"^. For all curves we find a small non- 
zero intercept whose value varies with vq. This introduces an extra free parameter in the linear 
fit corresponding to the non-zero intercept. The curve with the smallest non-zero intercept has 
ro = 2.5r] and gives a value of g — 0.47 with an error of the order of approximately 10% depending 
on the time range (here taken to be 15t^ < t < 75tjJ. This value of g is smaller than that found by 
Yeung and Borgas 0| and Ishihara and Kaneda , though still of the same order of magnitude, 
but agrees well with that of Ott and Mann fT3l| and Boffetta and Sokolov ■ 

In order to make a more complete analysis of Richardson's model, we compute the pdf of the 
separation distance. The Richardson pdf relies on two phenomenological assumptions: the first is 
that the eddy-diffusivity is self-similar, the second is that the velocity field is short-time correlated. 
However, it is known that anomalous corrections to the K41 scalings exist (see e.g. and these 
are likely to complicate the situation. 

In Fig. 121 we compare the separation pdf for the smallest initial separation, ro — 1.2?7, calculated 
from the DNS data, with that predicted by Richardson, namely For small times (up to 

t ~ 40t^), we observed a rapid change in shape with the pdf showing a pronounced tail, which 
indicates that while most pairs are still close together some have moved very far apart (not shown). 
At these times the curves do not rescale indicating that the early stages of the separation process 
is very intermittent. Here, the physics of the dissipative range still exerts an influence on the 
separation process and so we would not expect agreement with the Richardson pdf. Only for times 
in the range 40 — 70r^ do we find reasonable agreement with the Richardson pdf. We note that 
while at i ^ 40r^ we find good agreement for the tail but a large mismatch for values close to the 
peak, at t ~ 70t^ the pdf is almost undistinguishable from At large times the particles are 
moving more or less independently and so the pdf of will be a chi squared distribution with 
three degrees of freedom (not shown). 

A more detailed analysis of the separation pdf can be made by considering the separation skewness, 
Sr{t) = ((r-(t) -^)3)/((T2(t))3/2, and the kurtosis, ICr{t) = {{r{t) -Hi)f)/{al{t)f, where r is 
the mean separation distance and Ur is the root mean square separation distance. These are shown 
in Fig. Ofor tq = 1.2r] and clearly show the intermittent nature of the separation process at small 
times. The Richardson pdf, of course, predicts constant values for the skewness and kurtosis, 
namely 1.7 and 7.81 respectively and which are not reached until approximately t ~ 35r^. This 
time is within the inertial subrange and we may have expected the skewness and kurtosis to level 
off before decreasing to their large time values (0.49 and 3.1 respectively). That this is clearly 
not the case suggests that either contamination of the inertial range due to the dissipative and 
integral scales prevents us from having a region of constant skewness and kurtosis, or points to 
shortcomings in the Richardson model. 

These results put the difficulties of calculating Richardson's constant (described above) into con- 
text. A perfect collapse of curves in the pdf would have implied self-similarity but its absence is not 
necessarily an indication of the failure of the Richardson model; as we have already discussed, each 
end of the inertial range is affected by, respectively, dissipation range and integral scale effects. In 
section Fill BI we show how these problems may be overcome by measuring statistics at fixed scales. 

We conclude this section by considering the correlation function R{t, r) = {r{t)r{t + r)) of the 
separation distance for travel times within the inertial subrange. This quantity, which probes 
two different times along the separation process, is influenced by the temporal properties of the 
turbulent energy cascade sustaining the separation growth. In the inset of Fig. 21 we plot R{t,T) 
for — t < r < at different travel times t for pairs with initial separation rg = 1.2ry. In agreement 
with JuUien et al. p^ . we flnd that R{t,T) broadens with increasing travel time indicating that 
the pairs decorrelate more slowly at larger travel times, a consequence of the fact that larger and 
larger eddies have slower and slower dynamics. In the body of the flgure we plot the same data 
versus r/t. Dimensional analysis shows that R{t,T)/ R{t,0) — /{r/t). We see that all the curves 
collapse onto a single curve supporting the dimensional prediction. 
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B. Fixed scale statistics 



To disentangle the effects of different scales, an alternative approach, based on exit time statistics, 
has been proposed This consists of fixing a set of thresholds, r„ = p"ro, where p > 1 and 

n = 1,2,3, .. . and then calculating the time, T, taken for the pair separation to change from r„ 
to Vn+i- By averaging over the particle pairs, we obtain the mean exit time, {Tp{rn)), or mean 
doubling time if p = 2. Formally we are calculating the first passage time. The advantage of this 
approach is that all pairs are sampled at the same scales and that finite Reynolds number effects 
are less important |23. For particle pairs with initial condition p(r,0) = p'^S{r ~ rn/p)/4'Kr'^, a 
perfectly reflecting boundary condition at r = 0, and an absorbing boundary condition at r = r„, 
the pdf of the exit time, T, is defined to be 



'Pp,rAT) = -^ I P(r,r)dr, (4) 



from which we get 



7'p,.„(T) = -4^fcoei/Vr3 



dp 
dr 



(5) 



on making use of J^l- Following Boffetta and Sokolov we can derive a solution of the 3D 
diffusion equation in terms of an eigenfunction decomposition. This gives us 

oc 

= ^Qexp(-A?t)r'/V7/2(3A,0, (6) 

where ^ = {kQe^/^)~^/'^r^/^ , Jijiis^) is a Bessel function of the first kind. A; = 
1/3 (fco£^^'^)^^^?'~^^'^j7/2,i where ^7/2,4 are the zeros of J-!/2{x) and Ci are constants. It then follows 
from (jSJ that the large time asymptotic form of the exit time pdf is given by 

where k « 2.72 is a numerical constant derived from the leading zero of the Bessel function 
described above. 

Using Richardson's diffusion equation (Q, the mean exit time can be shown to be 

1 2/3 _ 1 2/3 

In the body of Fig. we plot {Tp{rn)) for a range of initial separations. It is immediately clear 
that there is no dependence on the initial separation in contrast to the mean square separation 
calculated as a function of time (see Fig. Moreover, we see a clear inertial scaling region in 
which the mean exit time grows almost like r^/^. 

Equation provides us with a method for calculating the Richardson constant (since fco is related 
to g): 

^ Ui {pV^-lf 
^ 81 p2 e(Tp(r))3 ' ^ ^ 

In the inset of Figure|Sl we plot the expression for the Richardson constant versus r, for various 
initial conditions. We see that a collapse of curves is beginning to form for all initial separations. 
We estimate the value of g to be approximately 0.50 ± 0.05, which agrees with the value computed 
above and with previous estimates of g [T^ . Il^ l29l | . This method has the advantage of relative 
insensitivity to the initial separation and avoids the problem of the non-zero intercept discussed in 
nil Al Of course, the present calculation of g assumes the validity of Richardson's model. We find 
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that g does not change significantly for p G [1.15, 2]. It is also worth noticing that the estimate of 
g is not very sensitive to the Reynolds number (see inset of Fig. |SJ). 



The exit time pdf, Vp.r^iT), is shown in Fig. for tq = 1.2r] and clearly shows the exponential 
nature of the exit time pdf at large times. Moreover, it shows a clear collapse of curves for 
intermediate and large exit times, Tp(r„) > (Tp(r„)), indicating that the exit time statistics in 
this range are self-similar. We note here that this collapse deteriorates with increasing p. We 
have shown that by focusing on properties at fixed scales a good agreement with the theoretical 
prediction © can be achieved. Free of the effects of the infrared and ultraviolet cut-offs, the 
Richardson diffusion model appears to work well for the inertial range of scales. For small exit 
times, Tp{r„) <^ {Tp{r„)), on the other hand, we do not find a complete collapse of curves at 
different thresholds, indicating that pairs with a fast separation are likely to exhibit intermittency 
(see inset of Fig. 01. 



The higher order moments of T are dominated by those pairs which separate slowly. Conversely, 
the moments of the inverse exit times, {[l/Tp{r)]P), are dominated by those pairs which separate 
rapidly and correspond to positive moments of the separation. Kolmogorov scaling based on 
dimensional analysis then leads to 
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eP/V-2p/3. (10) 



Assuming that a reasonable estimate of the exit time is given by T(r) ~ r/ur, where Ur is the 
relative velocity at scale r, intermittency corrections can be quantified in terms of the multifractal 
formalism |3Q||: 



Tp{r)J I Tl \Lo 



^ ^ Ce{p)-p 



(11) 



where Ce{p) are the scaling exponents of the Eulerian velocity structure functions as predicted 
by the multifractal formalism. In Fig. |7|we plot {{l/Tp{r))'''y^P scaled by the Kolmogorov (|10|l 
and intermittent scaling exponents l(TT|l . respectively. The Ce{p) are calculated using the She- 
Leveque formula [sj- As already remarked at lower Reynolds numbers by Boffetta and Sokolov 
|14| , there is a small but clear improvement in the scaling of the inverse exit times when scaled by 
the multifractal predictions. 



Before concluding this section we note that the exit time statistics can be used to measure the 
largest Lyapunov exponent in the flow. This is because for small thresholds, r„, the mean exit 
time probes the exponential growth of the separation distances. The exact relation between the 
'finite size Lyapunov exponent' and the mean exit time is 

In Fig. IHlwe show the right hand side of H12|l for three different Reynolds numbers (two from this 
numerical simulation, see Table 1) and one from a previous DNS study at different thresholds, 
r„. The usual Lyapunov exponent is recovered from the saturation value in the limit of small r„. 
As may be seen in the figure, the data show a clear proportionality between the Kolmogorov time, 
Tri, and the Lyapunov exponents. A, for all available Reynolds numbers. Thus, we get 

XTri - 0.115 ±0.005. 

This value is comparable with the one found by Girimaji and Pope 
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IV. RELATIVE VELOCITY STATISTICS 



A. Fixed time statistics 



We now consider the statistics of the relative velocity of the particle pairs during the separation 
process and which we denote as Ur{t) = u(^^(t) — u*^^)(t). The relative velocity statistics are 
of interest because they provide information on the rate of separation of the particle pairs. We 
consider the statistics of the relative velocity projected in the direction of the separation vector, the 
'longitudinal' component, and the projection of the relative velocity orthogonal to the separation, 
the 'transverse' component. The former is given by 

dr 

^ll=-=u..r, 

where f — r/r. The transverse component of the relative velocity is given by 

There are, of course, two transverse components of the relative velocity but since the turbulence is 
isotropic it suffices to consider only one. We comment here that the relative magnitudes of (|ur|), 
and (|uj_|) and the ahgnment properties of u^, r(i) and r(0) have been discussed extensively 
by Yeung and Borgas 0. Here, we state simply that our data give similar results and concentrate 
on the pdfs of the velocity components and their properties. 

In Fig. Owe plot the pdf of the longitudinal component of the relative velocity, for tq — 1.2r/. 

The pdf is negatively skewed at t = (not shown), corresponding to the Eulerian distribution, 
but as t increases, it quickly becomes positively skewed indicating that pairs with small initial 
separation are more likely to be diverging than converging. This skewness then decreases and the 
pdf tends towards a Gaussian distribution for travel times of order Tj, . The pdf of one component 
of u_L for the same initial separation is shown in Fig. 1101 Unlike the pdf of it is symmetric 
about the mean. Thus, negative velocities are as common as positive velocities indicating that 
there is no preferred direction of rotation as may be expected in isotropic and parity invariant 
turbulence. We note here that for both longitudinal and transverse pdfs we do not see a complete 
collapse of curves for times in the range t e [10, 70]t^. 

We consider the pdfs of the relative velocity components in more detail by analysing their skewness 
Su(t) and kurtosis /C„(i). These are shown in Fig. ^2 for rp — 1.27]. At t = the Lagrangian 
statistics (not shown) are identical to the Eulerian statistics. This is reflected in the negative 
skewness of U|| which is close to —0.55, the value commonly observed for Eulerian velocity structure 
functions at moderate to high Reynolds numbers jlH. We also note that at early times (up to 
t ~ lOr) the maximum values of the skewness and kurtosis of the velocity statistics are higher than 
the corresponding maxima of the separation statistics. We conclude this section by measuring the 
correlation of the relative velocity along the particle pair trajectories. In the inset of Fig. 1121 we 
plot D{t,T) = (M||(t)u||(t + r)) for — i < t < for pairs with initial separation tq — 1.2r]. In 
agreement with Fig. 21 we find that D{t,T) broadens with increasing travel time confirming that 
the velocity decorrelates more slowly at larger travel times. In the body of the same figure we plot 
the same data rescaled versus r/t. We note here that the rescaling does not give as good a collapse 
as for the separation distance, a fact which suggests the presence of tiny anomalous fiuctuations 
in the characteristic times governing the decorrelation of eddies of different size. 



B. Fixed scale statistics 



Following the exit time method of section UlIBI we calculate the relative velocity at fixed scales in 
order to achieve 'uncontaminated' inertial range statistics and which we term the exit velocities. 



8 



We compute the value of the relative velocity components W||(r) and uj_(r) whenever a particle 
pair has a separation within a specified logarithmic shell of radius r = r„(l ±0.1), with r„ = p^rg. 
This differs from the method we used to calculate the exit times above as here we are calculating 
not just the velocity at the first passage but also at all subsequent passages. 

In Fig. ^] we plot the mean longitudinal exit velocity, (u||(r)), as a function of the absolute 
separation, r, for different rp. The lack of dependence on ro is immediately apparent. We also 
see a clear separation between dissipative and inertial range scales, with growing linearly for 
small r, corresponding to the dissipation range and then growing close to r^/"^ for larger r which 
corresponds to the inertial range. Similar behaviour was observed for |uj_(r)|. We note here that 
the statistics showed little qualitative change for different values of p. It is tempting to think that 
because of the lack of dependence on the initial separation and the qualitative agreement with 
the K41 scalings for Eulerian statistics, the exit velocities resemble the Eulerian velocity structure 
functions even if they are not strictly Eulerian. However, the very existence of a non-vanishing 
mean longitudinal relative velocity tells us that the two sets of statistics cannot be identical. We 
consider these differences in more detail by analysing the pdfs of the exit velocities. 

In Fig. El we plot the pdf of the longitudinal (exit) velocity, P{u^\), and in the inset we plot 
M|l'P(u||). In contrast to the Eulerian pdf, the pdf is slightly positively skewed initially but as 
the separation threshold increases, it becomes more symmetric and tends towards a Gaussian 
distribution. However, unlike the pdf of M||(t) (see Fig. O we do not see an initial rapid increase 
in the positive tail of the pdf. At fixed separations there is not the spread of contributions to the 
velocity as there is at fixed times. The pdf of one component of the transverse relative exit velocity, 
V{u±z), is shown in Fig. 1151 for the same initial separation tq = 1.2?^. As may be expected, this 
pdf is initially symmetric and remains so for increasing separation threshold. We note that there 
is no qualitative change in the behaviour of both pdfs for different values of p. The absence of a 
complete collapse of curves at different thresholds in both pdfs indicates that the exit velocities 
are intermittent. 

We examine the intermittency of the exit velocities by considering their second and fourth or- 
der moments. Since the exit velocity statistics resemble Eulerian velocity statistics, we use the 
multifractal formalism for Eulerian velocity structure functions to quantify the intermittency cor- 
rections. In this way, we have a reasonable estimate to compare with. The Eulerian structure 
functions scale according to (Srv) (r/Lo)^^^^^ (see e.g. jl^). The She-Leveque formula for the 
scaling exponents, Ce{p), gives C-e(2) = 0.7 and Ce(4) = 1.28. Fig. [TCI shows that the second or- 
der moment of both the longitudinal and transverse relative velocity components scales according 
to the multifractal prediction. For the fourth order moment we find that the transverse compo- 
nent scales well with the multifractal prediction but the longitudinal component shows a small 
discrepancy. 



V. CONCLUSION 



We have considered the separation process of particle pairs in homogeneous isotropic turbulence 
in considerable detail. In addition to presenting such classical but important statistics, as the pdf 
of the separation distance and its second order moment (which gives us Richardson's constant), 
we also considered higher moments of the separation. Here, we found that dissipation range and 
integral scale effects made a quantitative assessment of intermittency corrections in the inertial 
subrange difficult. In order to get a clearer separation between the dissipation range, inertial 
subrange and integral scales we computed the statistics at fixed separations, the exit time statistics. 
This provided us with an alternative method for calculating Richardson's constant whose value 
agreed well with the 'direct' method. Moreover, these statistics allowed us to estimate intermittency 
corrections quantitatively in terms of the multifractal formalism. In agreement with Boffetta and 
Sokolov we found a small but noticeable deviation from self-similar behaviour for those pairs 
that separate rapidly. Hydrodynamic turbulence is most likely to be responsible for this rapid 
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separation of a small number of pairs; support for this comes from the improved scaling that the 
multifractal model - an inertial subrange model - provides for those particles that separate rapidly. 

We also calculated the longitudinal and transverse components of the relative velocity as a function 
of both travel time and the separation, that is, the exit velocities. Analogous to the separation 
statistics, we found that dissipation range and integral scale effects made a quantitative assessment 
of intermittcncy corrections difficult in the inertial subrange for the Lagrangian statistics computed 
as a function of time. However, the fixed scale approach showed a clear separation of scales with 
the velocities scaling like the Eulerian velocity structure functions (for moments up to order four). 
These statistics allowed us to quantify intermittency corrections in terms of the multifractal model 
for Eulerian velocity structure functions. However, we also noted a small but significant difference 
with the true Eulerian statistics for the case of the longitudinal exit velocity. 
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Rx Urms e V T] Lq Te Tl Tr, T dx Np 

183 1.5(1) 0.88(8) 0.00205 0.01 3.14 2.1 1.3 0.048 5 0.012 512^ 0.96-10'' 
284 1.7(1) 0.81(8) 0.00088 0.005 3.14 1.8 1.2 0.033 4.4 0.006 1024^ 1.92 •10'' 

TABLE I: Parameters of the numerical simulations: Taylor scale Reynolds number R\, root-mean-square 
velocity Urma, mean energy dissipation e, viscosity v, Kolmogorov length scale 77 = (i^^/e)^^*, integral scale 
Lo, large-eddy turnover time Te = Lo/urms, Lagrangian velocity autocorrelation time Tl, Kolmogorov 
time scale — {u/e)^^^, total integration time T, grid spacing dx, resolution TV'' and the number of 
Lagrangian tracers Np. 
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FIG. 1: The evolution of {r{t)^)/rj^ vs f/r,, for the initial separations ro = I.277, ro = 2.577, ro = 9. 877, and 
ro = 19.6?;. The straight line is proportional to i!^ . Inset: {r{t)^)/et'^ for the same four initial separations 
starting from t/r^ ~ 15. 
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FIG. 2: Comparison of the Richardson pdf with the DNS data. The curves refer to data for ro — I.277 at 
t — 5.2r^ (sohd line), t — 7r,, (long dashed line), t = 14r,, (short dashed hue), t = 42r^ (dotted line) and 
t = 70Tri (dot-dashed line). The thick solid line is the Richardson pdf ij^Jl. 
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FIG. 3: The separation kurtosis for the smaUest initial separation ro — 1.2rj. Also shown at some times 
are the error bars calculated from the minimum and maximum values of the five sub-ensembles. Inset: the 
separation skewness for the same initial separation together with the error bars at the same times. The 
horizontal lines are the appropriate values derived from the Richardson pdf and the chi squared distribution 
with three degrees of freedom. These values are 1.7 and 0.49, respectively, for the skewness and 7.81 and 
3.1, respectively, for the kurtosis. 
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FIG. 4: The normalized correlation function R{t,T) / R{t,0) versus r/t for ro = 1.2r; at different travel 
times. Inset: the same correlation functions now plotted versus t/r^. Curves (from left to right) refer to 
travel times t — 77t^, t — 63t^, t = 49r^, t = 35t,, and t = 21t^. 




FIG. 5: The mean exit time for the initial separations ro — 1.2rj (thin continuous line), ro = 2.5j7 (long- 
dashed line), ro — 9.8r] (short-dashed line) and ro = 19.6?; (dotted line) with p — 1.25. The straight line 
is proportional to r^^^. In the inset we show Richardson's constant, g, versus r/rj as given by (|5J for the 
same initial separations at R\ = 284. To evaluate the variability of g with the Reynolds number, we also 
plot a curve (thick continuous line) for the initial separation ro — 1.2j7 at Rx — 183. 
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FIG. 6: The log-linear plot of the exit time pdf for ro = l-2rj with p = 1.25 at r = 21.8?7 (dashed line), 
r = 83.3?) (dotted line), and r = 130.1?? (dot-dashed line). The solid line is the large time prediction Q. 
Inset: a lin-lin plot of the same figure showing more detail. 
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FIG. 7: The inverse exit time moments, {[1/Tp{r)]^)^''^ , for p = 1, . . . , 4 compensated with the Kolmogorov 
scalings (solid lines) and the multifractal predictions (dashed lines) for the initial separation ro — 1.21) and 
for p = 1.25. 
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FIG. 9; The pdf of W|| for ro = I.277 at the following travel times: (from outer to inner curve) t — 5.2r^, 
t — iTr), t = 14r,,, t — 42r,, and t = 70r^. The thick solid line is a Gaussian distribution. 





FIG. 11: The kurtosis of U|| (continuous line) and of u±_ (dashed line) for the smallest initial separations 
ro = 1.2»7. Inset: the skewness of M|| for the same initial separation. The horizontal lines are the Gaussian 
values for the kurtosis and skewness. 
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FIG. 12: The normalised correlation function D{t,T)/ D{t,0) versus r/t and t/r,, (inset) for pairs of initial 
separation ro = 1.2r). Curves (from left to right) refer to travel times t = 77t,,, t = 63r,,, t = 49r,,, t = 35t,, 



and t = 21t,,. 
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FIG. 13: The mean longitudinal exit velocity, as a function of the separation and scaled by w,,'r?/r,,, 

for pairs of initial separations ro = l-2?7, ro = 2.5r], ro = 9.8ij and ro = 19.6?; with p = 1.25. The thick 
solid line is proportional to r and the dashed line is proportional to r^^*. 
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FIG. 14: The log-hn plot of the exit velocity pdf 'P(u|[) calculated as a function of r for pairs with ro = 1.2r7 
and p = 1.25. The curves are for the following thresholds r = 5.72?7, r = 21.8?7, r = 83.3?7 and 130. Ir? (from 
outer to inner curve). The thick solid line is a Gaussian distribution. Inset: a lin-lin plot of W|'|7-'(w||). 
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FIG. 15: The log-lin plot of the exit velocity pdf Viu^^) for ro = 1.2r; with p = 1.25, at the same thresholds 
as Fig. 1141 The thick solid line is a Gaussian distribution. 




FIG. 16: The second and fourth order exit velocities for the longitudinal uy (x symbol) and transverse 
u_Lz (+ symbol) relative velocities for ro = 1.2rj with p = 1.25. The dot-dashed line are the multifractal 
prediction for the second and fourth order moments. 



